Theory of Vibration of the Larynx 1 

By R. L. WEGEL 

The vibration in the larynx is caused by an automatic modulation by 
the vocal cords of the air stream from the lungs. Analytically the mechan- 
ism is the same, and physically, closely analogous to that of the vacuum 
tube oscillator. It depends principally on the resonance of the vocal 
cords, the modulation of air friction in the glottis by their motion and 
the attraction due to constriction of the air stream between them. When 
these foices exist in certain relative proportions and phases, sustained 
oscillation as in singing takes place. The whole mechanism may be rep- 
resented analytically by force equations, from which conditions for accre- 
tion or subsidence of the vibration or for sustained oscillation may be 
easily deduced. These equations also show the analogy with other types 
of oscillating systems. 

IT is customary in treating the theory of the voice to assume the 
glottis or space between the vocal cords to be a source of a steady 
stream of air with superimposed periodic impulses caused by the 
vibration of the vocal cords. The harmonic content of these impulses 
is modified by the "resonating" vocal cavities before being radiated 
into free air. It is the nature of this modification which receives most 
attention. The mechanism by which the vibration of the vocal cords 
is maintained has not been carefully studied. 

The vocal cords are maintained in a state of sustained vibration 
by the proper balance between the various mechanical constants of 
the complete system, which thus act as a transformer of a part of 
the non-vibratory power derived from the air stream from the lungs 
into the vibratory power resulting in sound. It is a simple theory 
of this mechanism which is considered here. 

The method used is to obtain the force equations, which describe 
the vibrations of the complete mechanical system, by means of the 
Lagrange equations, from expressions of the total instantaneous 
kinetic and potential energies, the instantaneous forces acting and 
rate of dissipation of energy. The resulting simultaneous equations 
relating to the displacements and velocities of the various parts are 
then studied to find the frequencies of free vibration and the relations 
which must obtain between the various mechanical parameters of 
the system in order that one of these frequencies be sustained. The 
method is an application of the theory of H. W. Nichols, published in 
Physical Review, August, 1917. 

The theory is reduced to easily workable form by the introduction 
of simplifying approximations which will be described in the progress 

1 Presented before Acoustical Society of America, May 11, 1929. 
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of the discussion. The principal one of these is the neglecting of 
all reactions of second or higher order, thus leaving a set of linear 
differential equations. 

Structure of the Vocal Tract 

The vocal tract consists of three principal parts, the lungs and 
associated respiratory muscles for maintaining a flow of air, the 



FALSE VOCAL 
CORDS 



EPIGLOTTIS 




GLOTTIS 
VOCAL CORDS 



•-TRACHEA 



Kig. 1 — Anterior-Posterior Section of the Larynx. 

larynx (see Fig. 1) for producing the periodic modulation and the 
upper vocal cavities, pharynx, mouth and nose for varying the rela- 
tive harmonic content of sound originating in the larynx. 
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The capacity of the lungs in an adult man is capable of being varied 
from about two to five liters. The average in quiet breathing is 
about 2.6 liters. The average expiration of air in quiet breathing is 
about .5 liter. The rate of expiration of air in medium loud singing 
varies from 40 to 200 cm. 3 /sec, the lower values obtaining for trained 
singers. 

The larynx (see Fig. 1) consists of an irregularly shaped cartilaginous 
box at the top end of a tube, the trachea, about 12 cm. long by 2 cm. 
in diameter, leading from the lungs. The larynx contains the vocal 
cords, a pair of fibrous lips which in vibrating vary the width of the 
slit called the glottis, between them. The length of the glottis in 
the adult male averages about 1.8 cm. and in the female 1.2 cm. 
The width of the glottis varies widely with differing sounds. A few 
tenths of a millimeter may be considered representative. The tension 
and separation of the vocal cords are controlled by muscles. 

The principal upper vocal cavities are the pharynx, a space just 
over the larynx, the mouth and the nasal cavities. The first and 
second may be varied in size and shape at will, but the effect on the 
last is controlled only by varying the communicating aperture be- 
tween it and the pharynx. 

Equations of Motion of the Larynx 

Fig. 2 shows a cross-section of a model which illustrates the 
essential details of the larynx in so far as it is necessary for this treat- 
ment. kS'o represents the area of the opening to the trachea. The 
vocal cords are represented by elastically hinged members of com- 
bined effective area So. By effective area is meant the area of aper- 
ture which displaces the same volume of air as the vocal cords when 
it moves the distance 52 of the tips of the cords. This area is less 
than that of the vocal cords. The tips of the vocal cords are separated 
to form a gap, the glottis, of area Si. A positive or up and outward 
displacement g 2 of the vocal cords increases Si. It will be assumed 
that the air is not appreciably compressed in the neighborhood of 
the glottis, that is, any tendency to compression is relieved by flow 
into the trachea or pharynx. 

The pressure in the lungs forces a steady current of air through 
the glottis. Let the velocity in the trachea of this steady flow be Iq 
and in the glottis I\, Small vibrations of the vocal cords superimpose 
additional small velocities, i and i\, in the trachea and glottis re- 
spectively. If the instantaneous velocity of the vocal cords be h 
and it be assumed that they are constrained to move in synchronism 

(Jo + i )S = (7i + h)Si + i 2 S 2 . (1) 
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The above material is a description of a simple model of two degrees 
of freedom which simulates the principal characteristics from the 
standpoint of performance of the more complex larynx which has 
many degrees of freedom. It is this idealized model which will be 
considered in the subsequent treatment. Such points of performance 
of the actual larynx which may be due to the action of ignored and 




Fig. 2 — Schematic Larynx Model. 

presumably subsidiary modes of motion will, of course, not be pre- 
dicted by the theory. These are assumed to be of minor importance. 
The possible independence of motion of the two vocal cords will be 
considered later, however. 

The contraction of the air stream at the glottis introduces a rela- 
tively large concentrated kinetic energy in the air stream at this point 
similar to that at the mouth of a Helmholtz resonator. The inertia 
of a small plug of air between the vocal cords may then to a first 
approximation be treated as a mass L\. A concentration of frictional 
resistance also occurs at this point due to viscosity and to turbulence. 
A positive displacement q* (outward) of the vocal cords causes an 
increase in the mass of the plug of air in the glottis and a change in 
the effective resistance, R, encountered by it. The inertia L\ and 
resistance R of the glottis are therefore both functions of q», the dis- 
placement of the vocal cords from a mean position, and of the width 
of the glottis. If further Q 2 represent the average displacement of 
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the vocal cords from an appropriately chosen reference position, L 2 
their inertia coefficient and K» their effective stiffness, all measured 
at the tips, the total kinetic energy, T, and potential energy, V, of 
the larynx are 

T - §Za'(Ii + i\Y + |L2« 2 2 , (2) 

V = \K,{Q«. + g-i) 2 . (3) 

The Lagrange equation of forces for the nth coordinate of any 
system is 

dt di n dq n dq n , 

in which F n is a reaction due to friction. The force equations for 
the glottis and vocal cords therefore become 



Eo = Fl+ Ttdir 



WVx + *0* 



(5) 



0-Z^+* + *(ft + rt-<^^. (6) 

Nature of the "Constants" of the System 

It is quite safe to conclude that none of the coefficients (inertia, 
dissipation and stiffness) of the larynx are sensibly constant over the 
range of operation of the coordinates. Direct measurements are 
evidently impossible. It is conceivable that they may be arrived at 
indirectly by means of a comparison of experimental data, especially 
taken for the purpose, on voice curves and the results of dynamic 
analysis of the kind described here. The problem may also be studied 
by means of models. In order to solve equations 5 and 6 it is, how- 
ever, necessary to evaluate the space and velocity derivatives. 

A few simple experiments were performed on models for the sole 
purpose of determining the qualitative nature of variation of resistance 
of the glottis with displacement of the vocal cords. A diagram of 
the model used in the measurements is shown in Fig. 3. This con- 
sists of a brass tube, a, } / j t " in diameter, beveled off on the top at an 
angle of 45° with the axis, and two } $" brass plates, b, fitted on these 
beveled surfaces so as to leave a slit, S, which was made adjustable 
in width. A cross-section of this model is shown in c. The bottom 
of the tube was attached to a large air chamber in which the pressure 
and velocity of air flow could be regulated and measured. 

Three shapes of "glottis" were measured. The first had square 
corners, as shown on Fig. 3c. The second, 3d, was the same as 3c, 



212 



BELL SYSTEM TECHNICAL JOURNAL 



except that the corners of the lips were rounded. The third, Fig. 3e, 
had square corners as before, but the slit was about .1 mm. wider 
in the middle than at the ends. 




Fig. 3— Glottis Models. 

The resistance R is given as the ratio of the product of pressure 
and slit area to the linear velocity of flow. Measurements were made 
in each case through a range of pressures such as to give fluxes through 
the slit through a range of 50 to 200 cm. 3 /sec. (Stanley and Sheldon 
values, see Sci. Am., Dec. 1924) and through a range of slit width W 
of .01 to .10 cm. The data can be represented approximately in 
this range for the three slits by the following formula': 

R = 3.6 PW 1 - 7 X 10-°, 
R = 6.1 PW 2 - 2 X 10" 6 , 
R = 800 PW 2 - 3 X 10- G . 

In these expressions / is the velocity of flow of air through the slit. 
More careful data taken through a wider range of I and W would 
undoubtedly have given R in a power series. 

These formula? are taken to indicate that the resistance of the 
actual glottis increases faster than a linear function of / and W due 
to turbulence and may be represented as a single valued function of 
either displacement of the vocal cords g» (or glottis width) or of air 
velocity as expressed by a Taylor's series as follows: 
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In this expression i? is the resistance measured in the reference 
position at which point the derivatives are taken, where i\ and q* 
are zero. The experiment mentioned above determines the signs of 
the coefficients of q* and i\ as positive. If the flow were purely lami- 
nar, i.e. due to viscosity only, the first would be negative and the 
second zero. 

The Reaction F\ 

By definition, E = RqI\, where E is the force of the lung pressure 
on the glottis and 7i a corresponding linear velocity of flow of air. 
If a force Fi slightly greater than E act on the glottis and result in 
a velocity I = Ii -\- i it 

A, = R - < 8 > 

A combination of (7) and (8) constitutes an evaluation of F l for sub- 
stitution in the force equation (5). To a first order approximation 
then: 

doR \ . , T d R 



Fi = R h + {R + 7,^- ) k + 7, ^q,. (9) 

The coefficient of q 2 is dimensionally a stiffness and that of i x a re- 
sistance. In what follows they will be denoted by 

F, = Roh + Riii + K u q 9 .. (10) 

Glottis Mass (£,/) Reactions 

The kinetic energy of the air stream being proportional to the volume 
integral of the square of the velocity is largely concentrated in the 
glottis on account of the relatively high velocity at this point. On 
account of the irregularity in shape and turbulence in the stream it 
is impracticable to attempt an integration. If the velocity were so 
small that the turbulence were absent an approximate value of the 
air mass would be obtained by taking the mass of a cylinder of air 
having the length of the slit and a diameter equal to its width. This 
would make the mass L\ proportional to W 2 , or since the width is 
proportional to displacement of the vocal cords, to </ 2 2 - Owing, how- 
ever, to turbulence and other non-linearities, the mass is probably 
more nearly described as a tongue of air issuing from the glottis, 
the inertia L\ of which varies as some power function of the width 
and also of the velocity. 2 

- It has been found since experimentally (hat the mass reaction is very nearly 
that of a cylinder as described but reduced somewhat in diameter due to viscous 
or turbulent drag at the tips of the vocal cords. 
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It might be seen by carrying through an expression for this glottis 
mass involving a function of velocity similar to that for R of equa- 
tion (7) that only a quantitative change in effective mass would 
result in the final equations and that no new type of reaction would 
be introduced. This demonstration is not included here. In order 
to save space in this qualitative treatment it is ignored. For small 
displacements q» from a reference position at which the velocity of 
the air is I , the glottis inertia may be represented by the direct 
function: 

Ll Ll+ "^ ?2 + i a^ e2 " + etc - (11) 

in which the coefficient of q 2 is obviously positive. The second term 
of the second member of (5) may now be evaluated by performing 
the differentiations as indicated. Neglecting second and higher order 
terms and denoting dqo/dt by h the reaction in question becomes 

Lj£+h¥±H. (12) 

at dq 2 

The glottis mass of air, therefore, introduces two kinds of reactions: 
a simple inertia and a reaction proportional to the velocity of the vocal 
cords. For simplicity of notation (12) will be written 

L.f + %. (13) 

This completes the evaluation of the terms (5), the force equation 
of the glottis, which may now be written 

£o = Roll + Rxii + K u q, + M± + Gi,. (14) 

Force Equation of the Vocal Cords 
The force equation (6) of the vocal cords contains four terms. The 
first is the inertia reactance of the vibrating lips. The mass L 2 is 
the effective vibrating mass which, if multiplied by one-half the square 
of the velocity at the cord tip, gives the kinetic energy of their motion. 
If the distribution of the velocity in the vocal cords were known 
this might be found by integration. The second term F 2 in equa- 
tion (6) represents the internal dissipation and is assumed propor- 
tional to the small velocity i 2 . The third term is the elastic reaction 
which is proportional to displacement. 

The fourth term is a " gyrostatic " term. This term may be written 
as follows: 
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: -: '^+^ + etc -)- (15) 



dT i ,, . . ,„ I d L] 

t- = - J Ui + h) 2 -7— 
dq» \ dq« 



Again by neglecting second and higher order effects this reaction 
becomes 

±i doLi T doLj ■ I\~ dn~L\ t*£\ 

2 dg 2 3 go 2 dq-r 

It will be seen that the first term of this expression represents a 
static force tending, since it is negative, to draw the vocal cords 
together. This is the Bernoulli effect utilized in a venturi meter. 
This steady force is counterbalanced by an elastic reaction of the vocal 
cords with which it combines to determine an equilibrium position 
which obtains when the cords are not vibrating. This term may, 
therefore, be dropped from the final equations representing only 
superimposed motions. 

The coefficient of i\ is identical, except for a sign, with G of (13). 
It represents a force on the vocal cords due to a superimposed part 
of the Bernoulli effect caused by the small superimposed velocity i\ 
in the glottis. The coefficient of qi is dimensionally a stiffness. This 
apparent stiffness is due to the nature of the air flow and is inde- 
pendent of any elastic members. It is negative if the second differ- 
ential of glottis mass with respect to cord displacement is negative, 
positive when this coefficient is positive and vanishes when this 
coefficient is zero. It simply adds or subtracts in effect from the 
stiffness K 2 of the vocal cords. The first possibility is the more likely. 11 
These terms may then be written for simplicity 

- t- = - F - Gii- K x q«. (17) 

oqi 

Force Equations of the Larynx 

The force equations of the glottis and vocal cords with constants 
thus evaluated are 

E = Zi~ + Riii + Gk + RJ, + K u q 2 , (18) 

= Z 2 ^y + Riii + A' 2 go + K 2 Q 2 - F - Gii ~ K x q 2 . (19) 

As explained before, E = R Ii and F = K«Q«; so these cancel and 
are of no interest here. In the following it will be seen that the field 

3 This coefficient has since been found to be negative. 
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stiffness K x is included in K-> to simplify notation. This leaves (18) 
and (19) finally: 

= U ^ + *i»i + Gi 2 + K u q»_, (20) 

= L 2 ^ + Rih + K 2 q 2 - Gii. (21) 

It should be noted that these equations represent all first order 
internal reactions of the idealized model of the larynx. The series 
expansions have been carried out, to show to what approximations 
these equations hold. It should also be pointed out that the effects 
of mechanical hysteresis of the parts, which make the relative posi- 
tions of the parts dependent on the previous history of their motion, 
have not been considered. A consideration of hysteresis complicates 
the theory considerably and is ignored for the same reason and with 
the same justification and limitations that it is ignored in the ele- 
mentary treatment of electrical circuits containing coils with magnetic 
material and condensers with electrostatic hysteresis. 

External Reactions of the Trachea and Vocal Cavities on 

the Larynx 

So far the modifying effect of the trachea and lungs, as well as the 
upper vocal cavities, on the motion have not been considered. Before 
using the equations it is necessary to evaluate these reactions and 
add them in their proper places. 

Imagine a weightless piston fitted into the trachea just below the 
vocal cords such that the volume of air thus enclosed in the larynx 
is so small in comparison to that of the trachea and lungs that its 
compressibility may be neglected. If the vocal cords are held- rigid 
and the plug or piston of air in the glottis is forced inward, a reaction 
in addition to the resistance and inertia of the glottis will be encoun- 
tered due to the impeding effect of the trachea piston, which impe- 
dance is determined by the constants of the lower chambers. If a 
small force / act on the trachea, causing a small velocity, i , and 
we assume linearity of response / = Z io where Z Q is a constant which 
may, due to a positive inertia reactance or a stiffness, contributed 
by air compression in the lungs, involve either a time derivative or 
integral of displacement. For the present consider it to be a gener- 
alized impedance operator. Due to the relative incompressibility of 
the air in the larynx, the volume displaced by the trachea piston is 
ioSo = iiSi. Since also the instantaneous pressure inside the larynx 
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is constant on all its walls, including the surface of the trachea piston 
fo/Sa = fi/Si. We then have 

f x = zMk. (22) 

This reaction due to the trachea must be added to those of the glottis 
given in (20). In like manner if the effective area of the vocal cords 
is Si a reaction f-> must be added to their force equation 

f, = zMk. (23) 

'JO 

Due to the steady component of air flow there is a static component 
of pressure tending to force the cords outward. This is counter to 
the static Bernoulli term and again, if second order effects of small 
quantities be neglected, serves only to alter the equilibrium position 
and may therefore be disregarded here. 

When the glottis plug of air is displaced inward a force is exerted 
on the vocal cords tending to move them outward which is relieved 
to a certain extent by a yield of the trachea piston. This force on 
the vocal cords may be shown by reasoning similar to that above 
to be 

Z %=k. (24) 

Since this part of the system is linear, the reaction between glottis 
and vocal cords through this channel is reciprocal so a force is exerted 
on the glottis when the vocal cords are displaced of 



SiS, . 

So 2 



Z.=sr?i. (25) 



It will be noticed that Si is a variable because of the variation in 
width of the glottis while vibrating. The effect of this variation in 
these terms is obviously second order since i\ is small and will therefore 
be neglected. 

The reactions of the upper cavities might be similarly added, but 
they are apparently relatively small and since they are at present not 
quantitatively known, are disregarded in the general equations be- 
cause of the increased complexity. Generally, however, Zo may be 
thought of as representing the additive effects of both upper and 
lower chambers. 

The complete force equations of the voice for small vibrations, 
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taking into account all major external as well as internal reactions, 
may then be written : 

= Lx-jT + -Ri*i + G *2 + K u q 2 + -FTT- *i H rr^— *2- (26) 

Of Jo" "JO 



= Lo-T- + #2*2 + -^2<Z2 - G*l + -TrV^S 



"(M1O2 • 



(27) 



These equations may be put in a somewhat simpler form by virtue 
of the fact that they are linear differential equations with constant 
coefficients. In such a case the time differential may be replaced 
by an algebraic operator p such that i = pq, di/dt = p-q, where p is 
of the dimensions and nature of a frequency 

= ( fU + pRx + &•!§}) 2i +(& + PZ O ^0 + #„) &, (28) 



M | -/' G + Mo§# J 



+ f ^L 2 + £#•> + A% + pZo^ 
The determinant of this system is (calling pZ = Y ) 

(?L 1 +pR l + Y ^(pG+ Ko^ + Xu) 



7) = 



22- (29) 



(30) 



D - 



(31) 



This determinant represents the complete reactions of the larynx and 
the external effects of communicating air chambers. 

If the effects of the air chambers be disregarded the system is 
represented by placing Y - 0, giving the simple form 

(p'Li + pRi) (pG + K u ) 

- pG (p-L 2 + pRo + K 9 ) 

Nature of This System 

The voice system represented by determinant (3) is very closely 
analogous to other vibrators, such as the microphone oscillator or 
door buzzer and the vacuum tube oscillator. The literature on the 
latter subject is now so extensive that the pointing out of the analogy 
should make the method of solution for sustained oscillation, as in 
singing, or for subsidence or accretion of the oscillation, as in speak- 
ing, clear to any one familiar with it. 
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Fig. 4a is a schematic diagram of a three-element vacuum tube 
oscillator circuit known as the "tuned grid" circuit. This is one of 
many kinds. The transformer coupling between the plate and grid 
circuit is represented by an auto-transformer. Fig. 4b represents 
the same circuit schematically but with circuit elements only. In this 
Rn represents that part of the resistance of the coil which belongs to 
the grid circuit and any other associated dissipation, Z, 2 is the in- 
ductance of the coil as seen from the grid mesh and K«_ the reciprocal 
of the combined tuning capacity across the grid and that of the grid- 
filament. It is the electrical stiffness or elasticity of the grid mesh, 
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Fig. 4 — "Tuned Grid " Oscillator. 

in other words. Ri is a plate-filament resistance ("a. a") and L\ 
that part of the coil in the plate circuit. M is the mutual inductance 
of the transformer which is not part of the mesh impedance of either 
plate or grid. The element K u is the " uni-lateral mutual impedance " 
(G. A. Campbell, 1914) between the plate and grid meshes and is 
numerically equal to fj.K 2 , where n is the amplification constant of 
the tube. Other internal tube impedances are as usual neglected. 
The impedance determinant may be written directly from the circuit 
diagram Fig. 3b. 

(p'L, + pR t ) {pHI + K u ) 
i?M {p-U + pR l + K 2 ) 



D = 



(32) 



The quantities on the principal diagonal of this determinant, that is 
the first and last elements, are as usual in a circuit determinant the 
mesh impedances while the others are the mutuals. The principal 
features of the analogy may be seen by comparison of determinants 
(31) and (32). Except for the thus far undefined external or trachea 
impedance the mesh impedances are the same, from which it appears 
that the glottis is analogous with the plate-filament path in the vac- 
uum tube and the vocal cords with the grid-filament path. The air 
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velocity in the glottis I x corresponds to the plate current. In the 
vacuum tube this plate current is modulated by varying charge, q 2 , 
on the grid. In the larynx the glottis air velocity is modulated by 
varying displacement, q 2 , of the vocal cords. The charge on the 
plate (again neglecting internal capacities except the grid-filament) 
causes no effect on the grid mesh and in the larynx the position of 
any element of glottis air has no effect on the vocal cords. The uni- 
lateral mutual impedance, K,„ is the same in both. 

The analogy breaks down at the point where the " feed back " part 
of the mechanisms is compared. The " feed back" is the bilateral 
part of the mutual impedance between the two meshes. In the 
vacuum tube circuit this is p-M, the mutual of the transformer, while 
in the larynx it is pG, the " gyrostatic " mutual. The latter is a type 
of element which does not occur in electrical circuits, arising as it 
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Fig. 5 — Tuned Grid and Wind Reed Circuits. 



does from a variation of a mass or inductance with a displacement. 
Inductance, being a function purely of the geometry of a circuit, 
can only vary with mechanical displacement and not with electrical 
displacement or charge. The gyrostatic mutual is common in the 
mechanics of rotating bodies whence it derives its name. It is also 
the mutual in an electromagnetic telephone receiver or relay be- 
tween the electrical circuit and the armature or diaphragm. 

In order to fix the rather useful concept of the analogy in mind, 
Fig. 5 is added showing the schematic circuit of the vacuum tube 
(5a) and a circuit diagram (5b), which represents determinant (31) 
the characteristic formulation of the dynamics of the larynx. Fig. 5b 
is represented by the conventions of an electrical circuit, except for 
the element G for which a different convention is necessary. The 
one taken here is that of a resistance enclosed in a rectangle. From 
(31) it will be seen to be similar to a resistance in its association with 
frequency p but different from resistance in that it occurs non-sym- 
metrically in sign in the determinant. It does not involve dissipation. 
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Its occurrence here is the simplest possible for when there are appre- 
ciable concealed or ignored modes of motion it may have the form 
of a generalized impedance containing at least one element of resist- 
ance, but will always be non-symmetrical as a whole in sign in the 
determinant. 

The use of the circuit for representing the mechanical system is an 
extension of an old but recently popularized method of studying 
mechanical or electrical vibrating systems by the help of analogy, 
one with the other. The extension consists in the explicit representa- 
tion by diagram of the gyrostatic mutual which makes the deter- 
minant unsymmetrical in sign and of the unilateral mutual which 
makes the determinant unsymmetrical in magnitude. Fig. 6 is a 
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Fig. 6 — General Wind Reed Circuit. 

diagrammatic representation of the more general system of deter- 
minant (30). This includes the external F reactions as well as the 
internal. 

Having thus described the extended method of analogy the follow- 
ing study of the larynx with the help of the circuit diagram of its 
determinant should be clear. 

Sustained Vibration of the Simple Larynx 

In vibrating, the vocal cords do not receive excitation of the fre- 
quency at which they vibrate. The source of power is in the air 
stream 7\ which enters the equations in K„, the unilateral mutual 
impedance. Since this is treated as a constant circuit or dynamical 

element this air stream may be ignored as a drive and the resulting 
15 
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vibration considered as the free oscillation of the system. The de- 
terminant (31) (or 30) is then used to determine the free frequencies 
and decrements of the system. The method is as usual to solve for p 
in the equation 

D = 0. (33) 

To simplify the demonstration the simple larynx without the load of 
the air chambers will be considered. Taking D of (31) then and 
expanding: 

p*LiU + f(LJU + URi) + P n -{LJ<, + RyR> + G 1 ) 

+ p(R x K, + GK U ) = 0. (34) 

If this be divided by L\L» and the uncoupled decrements and natural 
frequency defined: 

(35) 
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= 0. (36) 

One of the roots of this equation is zero and another is negative 
real since all coefficients are positive. This root is therefore the decre- 
ment of a mode of non-vibratory motion. The remaining two roots 
may be real, imaginary or generally complex, of the form 

A ±./o). (37) 

If it is found that A = 0, then an oscillation once started will be 
sustained. If A be negative then any existing oscillation must sub- 
side or if A be found positive then an impulse will start an oscillation 
which of itself increases in amplitude to a point where its violence 
modifies the constants to such an extent as to make A vanish, leaving 
a sustained oscillation, or negative leaving the oscillation to subside 
to a lower amplitude or completely if sufficient permanent changes 
have been made. 

If now (36) be written 

Ap* + Bp- + Cp + D = (38) 

and the first root (37) be substituted, two equations result, one from 
the real and the other from the imaginary terms, as follows: 

^A(A 2 - 3co ? ) + 5(A 2 - a, 2 ) + CA + D = 0, (39) 

A(3A- - co 2 ) + B(2A) + C = 0. (40) 
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Now the condition for sustained oscillation is that A = and if the 
value of co when this obtains be co then 

, D o C .... 

= 7? a A' 

or the condition for sustained vibration in terms of the constants is 

AD = BC. (42) 

In addition to this if use be made of the fact that in an algebraic- 
equation such as (36) the coefficient of p~ is the negative sum of all 
the roots then this coefficient is the real root. Let this be A and 
then (36) may be written 

/> 3 + /> 2 A + pm 2 + Aocoo 2 = 0. (43) 

The coefficient of p in (36) is therefore the square of radial frequency 
at which sustained oscillation will take place and this is seen to be 
higher than the natural frequency cu 2 of the vocal cords, the difference 
being increased when the damping of either mesh is greater or when 
the coupling mutual G is greater. 

It might be noted in passing that (43) is the free oscillation equation 
for any system which may be represented by a cubic equation and 
is not confined to the simple larynx. Such an equation always results 
when there is only one kind of reactive element in one of the meshes. 
It holds also for the tuned grid circuit. 

The condition for sustained oscillation to be fulfilled for the con- 
stants may from (42) be reduced to: 

^-^(w- 1 )' (+4) 

It is rather difficult to place a simple physical interpretation on 
this formula. The qualitative import of it may however be seen by 
substituting the values of G and K„ from (13) and (10) : 



K\K- 



u = lw [ 



Ry Lx 






The first term in brackets is in the nature of a resistance modula- 
tion constant, a fractional change in glottis resistance per unit dis- 
placement of the cords, to be designated by r and the second term 
similarly a glottis mass modulation constant, /. The quantity LJ X 
is the momentum of the air in the glottis. This equation is then 

RyR, = {L x h)\r - /)/. (46) 
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Thus it appears that the resistance modulation must always be greater 
than the mass modulation and when the difference is small the air 
momentum must be increased to compensate. Owing to the physical 
limitation in accuracy of continuous maintenance of adjustment in 
the larynx, if a large momentum is depended upon to compensate for 
a small modulation difference, an unsteadiness or instability is likely 
to result. It is common experience that it is impossible to produce a 
sound with the voice with less than a certain minimum intensity. 
This corresponds, with the most favorable adjustment of the modu- 
lation constants which are physically possible, to a minimum momen- 
tum of air from the lungs which satisfies (46). It will be evident 
that this interpretation must not be taken too seriously quanti- 
tatively. 

Subsidence and Accretion of Vibration of the Simple Larynx 

Oscillograms made of the speaking voice show that, among other 
things, the amplitude of the oscillation and the pitch are in a con- 
tinuous state of change. This is also true in singing but not nearly 
to the same extent. It seems therefore that in singing the adjust- 
ment of the voice system for sustained oscillation as described in 
(44) above is of major importance, while in speaking conditions for 
variation are of most importance. 

The principle of the investigation of variation is simple enough 
but in all but the most elementary systems the algebra involved is 
impracticably awkward. If by solving (33) directly for the roots of p, 
it be found that A is positive, then any existing vibration will tend to 
increase while if A is negative, then vibration will tend to subside. 
The algebraic difficulties arise in the general solution but these are 
largely obviated by making the assumption, which is most likely 
usually fulfilled in practice, that the real parts of the roots may be 
treated as small quantities when compared with the imaginary parts. 
A common frequency for a man's voice is 150 cycles per second for 
which co is 1000 in round numbers. The decrement of a telephone 
receiver is ordinarily 100 to 200 in open air. The decrement of a 
tuning fork is represented by a fraction. Judging from variations 
in amplitude in an oscillogram (from which of course decrements 
may not be read directly) it would seem reasonable to assume that A 
is small compared with w - The study of variation thus becomes an 
investigation of small departures from a condition of sustained oscil- 
lation, the reference condition being that critical adjustment for which 
the roots of interest of (33) are pure imaginary. 
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Suppose in (38) that A — 1 ; then without loss in generality: 

p» +Bp*+Cp + D = 0. (47) 

In such an equation the roots are continuous functions of the co- 
efficients. The same is true of their derivatives except at the one 
point where transition occurs from pure real to complex values. The 
values of the roots of interest in this connection are in their complex 
region at the point where the real part of the root passes through a 
zero value. This is the point at which free oscillation of the oscil- 
lating mode occurs, the values of the roots of this mode being as shown 
before, ± juo. 

If it now be supposed that one cause produces small variations, 
directly or indirectly on each of the coefficients and that the magni- 
tude of this cause be x, then : 

( 3, + a , + O g + ,g + ,g + <g.0. (48) 

The problem then is to determine dp resulting from any assigned 
cause dx when p = jW From (43) we have at this point B = A , 

C = o>o" and D = Aoa>o 2 - 

.Ao\dp _ ,dB , . dC , dD 

o) 



2 "° = \ l -J-)£-- »<3I + ■'"« 57 + 37 • (49 > 



This is the frequency (complex) variation equation taken in the 
neighborhood of free oscillation. 

When any readjustment of the larynx takes place all of the "con- 
stants " entering the coefficients undergo change, in particular those 
of the glottis K,„ R u G. Suppose for simplicity that one only varies, 
then this variation dK lt , dR u or dG may be taken as the magnitude 
of the cause dx. In particular if K„ vary, 

dB = = dC and dD/dx = G/LiL 2 , 
1 _.A £ \_#_ G 

J io ) dK u 2m-LxU ' ( } 

If in addition A be small compared with o> , 

, GdK u i .Ao \ 

dp = 7 n T 1+./— • (51) 

Zcoo L\Lz \ coo / 

This shows that if a condition of sustained oscillation is departed 
from by slightly increasing K„, an increase in the amplitude of vibra- 
tion begins which is proportional to the logarithm, since (p -f- dp) 
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is the exponent, of the increment dK u and the frequency (imaginary- 
part) of vibration increases slightly in proportion. If K„ were the 
only varying element the vibration would continue indefinitely to 
increase. 

If on the other hand K„ be assumed constant, the variation being 
in Ri, then it may be similarly shown that 



2w " 






(52) 



whence it appears that a small increase in glottis resistance dR\ (or 
rfAi) introduces a subsidence of vibration but an increase in frequency 
of oscillation as before. A decrease — dR\ of course produces the 
opposite effect. 

If the change be in G, it turns out that 



tp. dG 



2wd 2 LiLo 



a icx ] : ( ^ + 2co G 



(53) 



Here it appears that an increase in the gyrostatic mutual, G, may 
introduce either a subsidence or an accretion in amplitude but like 
the others makes for an increase in frequency of oscillation. 

Variation in other elements produces similar conflicting tendencies 
not only in damping but in frequency. 

The physical picture to be drawn from this is that in speaking the 
voice modulates from one amplitude and frequency to another by 
proper relative variations in adjustments in its constants, being con- 
stantly in a state of changing subsidence or accretion. It would 
seem also that the principal cause of change in frequency is in the 
vocal cords and that of amplitude variation in the glottis. Speaking 
is, in this respect, a more intricate process than singing. 

Other Types of " Feed Back" 

The detailed study of the larynx has so far been limited to the 
assumption that the " feed back" is entirely gyrostatic. This is of 
course actually not the case. How much influence is exerted by the 
general Y is difficult to estimate. 

If the trachea were a long tube but still shorter than a quarter 
wave-length of sound at the frequency of oscillation and rather smaller 
in diameter, and substantially open at the end the mass of the air 
in it would then be appreciable and Y in (30) would be written p-L Q . 
If in addition the gyrostatic term were negligible the system would 
then be exactly analogous with the tuned grid system and (32) rather 
than (31) should be the subject of detailed study. 
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If on the other hand the lungs acted substantially as a solid walled 
chamber of comparatively small size, the elasticity of the contained 
air would be represented by taking K for F . The surface area in 
the lungs is very large compared with a regular chamber of equal 
volume so considerable dissipation must be encountered by vibration. 
If this were the most important reaction F should have been replaced 
by pR . 

Unquestionably all three types of reaction enter. A more general 
treatment to include them is plainly not a subject for a short paper. 
It is interesting however to note that in the dynamical system of 
brass horns these latter Y reactions exert controlling influences. In 
this case the lips of the player perform the same function as do the 
vocal cords of the voice while the external load, the horn, corresponds 
to the pharynx, the reaction of which is the same dynamically as 
the trachea. In this case the frequency of the horn is that of sustained 
oscillation and not that of the lips. The same is true of the wood- 
wind, in which case the reed or reeds replace the lips or vocal cords. 
In these cases Y is proportional inversely to the hyperbolic tangent 
of the frequency or may be approximately represented by the im- 
pedance of an anti-resonant element. 



